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Abstract 

We consider quantum systems with a chaotic classical limit that depend on an external 
parameter, and study correlations between the spectra at different parameter values. In par- 
ticular, we consider the parametric spectral form factor K(r,x) which depends on a scaled 
parameter difference x. For parameter variations that do not change the symmetry of the sys- 
tem we show by using semiclassical periodic orbit expansions that the small r expansion of the 
form factor agrees with Random Matrix Theory for systems with and without time reversal 
symmetry. 
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1 Introduction 



One of the characteristic features of quantum systems with underlying chaotic dynamics lies 
in statistical fluctuations of their spectra. If the energy levels are scaled such that their mean 
separation is one then the statistical distribution of the levels of individual quantum chaotic 
systems are found to be universal in the semiclassical limit h — > and to agree with those 
of eigenvalues of random matrices [1]. The appropriate ensembles of random matrices depend 
only on symmetries of the system. Systems with or without time reversal symmetry (TRS) are 
described by the Gaussian Orthogonal Ensemble (GOE) and the Gaussian Unitary Ensemble 
(GUE) (in the absence of half-integer spin and other symmetries, which we will assume in the 
following) . 

Universality can be observed, however, not only in the spectrum of an individual quantum 
system, but also in the way in which the spectrum changes due to an external perturbation. 
In the following we consider quantum systems that depend on an external parameter whose 
alteration does not change the symmetry of the system and for which the classical dynamics is 
chaotic for any parameter value. Correlations between the spectra of these systems at different 
parameter values are also found to be universal functions of the parameter difference provided 
that the parameter is scaled in an appropriate way [2, 3]. The universal correlation functions 
again agree with those of Random Matrix Theory (RMT) and have been derived for the GOE 
and the GUE. The early developments of parametric spectral correlations are reviewed in [4]. 

One main approach to understanding the connection between quantum chaos and RMT 
has been the application of semiclassical methods. One convenient quantity to characterize 
correlations within the spectrum of an individual quantum system is the spectral form factor 
K{t), the Fourier transform of the two-point correlation function of the density of states. 
Within the semiclassical approximation it is expressed by a double sum over the periodic orbits 
of the classical system. An evaluation of this double sum in the diagonal approximation, which 
pairs orbits with themselves or their time reverse, yields the first term in the small r expansion 
of K(t) in agreement with RMT [5, 6]. Higher order terms are due to pairs of correlated 
periodic orbits. The origin of these correlations was identified and the next-order term obtained 
for uniformly hyperbolic systems in [7, 8]. The contributions of all correlated periodic orbit 
pairs that contribute to the small r expansion of K(t) were evaluated and summed up by 
combinatorial methods in [9, 10, 11]. Similar methods have been applied since to derive off- 
diagonal terms, for example, for the conductance [12, 13], the shot-noise [14] and the GOE-GUE 
transition [15]. 

For parametric spectral correlation functions agreement with RMT has been obtained within 
the diagonal approximation [16, 17, 18]. In this article we go beyond the diagonal approximation 
and derive off-diagonal terms for the parametric spectral form factor. For systems without time 
reversal symmetry we derive all terms in the small r-expansion in a closed form. In the GOE 
case the method allows the calculation of arbitrarily many terms in the expansion. We compare 
the expansion up to r 7 with the result of RMT and find exact agreement. One main reason for 
the universal result is that in the semiclassical limit H — > the relevant quantum fluctuations 
are due to very small parameter variations on the classical scale. One assumption is that the 
parameter dependence is in some sense typical. Specifically, we assume that the derivatives 
of the actions of very long periodic orbits with respect to the parameter have a Gaussian 
distribution [17, 18]. This excludes, for example, rank-one perturbations for which off-diagonal 
terms were calculated in [19]. 

In section 2 we introduce the parametric spectral form factor and in section 3 we state results 
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of Random Matrix Theory for it. In section 4 we consider its semiclassical approximation and 
the diagonal approximation, while in section 5 we derive the off-diagonal terms. For systems 
without time reversal symmetry the expansion is summed in section 6, and section 7 contains 
our conclusions. 

While writing up our paper the preprint [20], which is closely related to our work, appeared 
on the archive. Nagao et al. investigate parametric correlations that depend on a magnetic field 
difference, and obtain the universal results for the GUE case and the GOE-GUE transition by 
periodic orbit expansions. Our work is complementary in that we treat arbitrary parameters 
and consider also the GOE case. 

2 The parametric spectral form factor 

One way to characterize fluctuations in quantum spectra is to consider correlation functions of 
the density of states. For parametric correlations the density of states depends on the energy 
E as well as on a parameter X, and in the semiclassical regime it can be written as the sum of 
a mean part and an oscillatory part 

d(E, X) = 5 ( E - E n(X)) « d(E, X) + d osc (E, X) , (1) 

n 

where E n (X), n = 1, 2, . . ., is the n-th energy level as a function of the parameter X. The 
mean density of states in an /-dimensional system is given by d(E,X) ~ Q(E, X)/(2irh)^ in 
the semiclassical limit h — > 0. Q(E, X) is the volume of the surface of constant energy in phase 
space at energy E and parameter X. 

In order to obtain a universal parametric spectral correlation function one has to perform 
two unfoldings, one in energy E and one in the parameter X of the system. A new energy 
parameter is defined by 

E = N{E,X), (2) 

where N(E,X) is the mean part of the spectral staircase N(E,X) = dE' d(E',X). In 

terms of the new energy E the density of states has a mean value of one. The spectral statistics 
are evaluated in the semiclassical limit in an interval AE that is classically small but contains 
a large number of energy levels, i.e. it satisfies E 3> AE 3> 1. 
A new parameter X is introduced by [2, 21] 

X = [ X dX> a(X') , a(X') = V >n(X') 2 > , (3) 

where v n (X) = dE n /dX are the level velocities and the average is performed over the levels 
in the interval AE. X is an arbitrary parameter value at which X = 0. In terms of the new 
parameter X the level velocities have a unit variance. 

We may then define the universal two-point correlation function by 

M, h *) = + + f ) t- (*-2*-f)) w , (4) 

where d(E, X) is the density of states of the unfolded spectrum, and the average is performed 
over the energy interval AE as well as over a parameter interval AX. The relation to the 
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original density of states is given by 



~ ~ ~ dN(E,X) dE d(E,X) 

rf(E ' X) " dE M~d(E^j' (5) 
and in the semiclassical limit we find that 

D , , ( d ° SC (E + Js + %,X + £) d~ ) 

RM Ke-xy • (6) 

Equation (6) has been obtained by linearizing the unfolding equations (2) and (3), because 
x and i] correspond to small changes on the classical scale, since d is of the order ft~f and 
a = cr(X) is of the order of ft~Cf +1 )/ 2 (see equation (25) later). The term x/2pa takes account 
of the change of the energy when X is changed while keeping E fixed 



dE 
P= dX 



ON/dX 



dN /dE 



In the following, we will consider the parametric spectral form factor which is obtained by 
a Fourier transform of the parametric two-point correlation function 

/oo 
R 2 (ri,x)e- 2mr > T drj . (8) 
-oo 

3 Results from Random Matrix Theory 

The parametric two-point correlation function R2(rj,x) has been derived in the context of 
disordered systems [2, 22] for the GUE and the GOE. For the GUE case it is given by 

-i pi poo 

B$ UE (r t ,x) = ±J d\J dAi cos(ttcu(A 1 - A))e-^ 2(Al2 - A2)/2 . (9) 
After performing the Fourier transform in (8) to obtain the parametric form factor we arrive at 

-i pi poo 

K gve (t, x) = ^J dA J dAi e"^ 2 ( Al2 - A2 )/ 2 [<S(Ai - A - 2r) + <S(Ai - A + 2r)] . (10) 

Because r is positive and Ai > A the second delta function does not contribute. From the first 
delta function we get the relation 2r = A x — A. In the case r < 1, which we consider in the 
following, the domain of integration for Ai is reduced to 1 < A x < 1 + 2r, and we obtain 

K^(r, x) = \ r T dA, e 2 ^^) = ^ff^ ^r T<1 (n) 
2 J 1 2n 2 x 2 r 

For comparison with the semiclassical expansion we expand the sinh-function and define B = 
2tt 2 x 2 /k where k — 1 and 2 for the GUE and GOE cases, respectively. 



' ^2n^4n+l 



tf- M= e-£|^. (12) 

n=0 v 1 
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The parametric correlation function for the GOE case is given by a triple integral 

= /_; cuf ^ f cu 2 cos(M a - a.a,)) ^:^-^^^ 

x e -^ 2 (2A 1 2 A 2 2 -A 2 -A 1 2 -A 2 2 +l)/4 _ ^ 



Evaluating the Fourier transform to obtain the parametric form factor results in a sum of two 
delta-functions 

K goe (t,x)= fdA TdAi /°°dA 2 - (1 ~ A2 9 )(A ~ AlA2) 2 2 ^ 

V ; 7-i A J i (2AAiA 2 - A 2 - Ai - A 2 + l) 2 

x e -. 2 x 2 (2A 1 2 A 2 2 -A 2 -A 1 2 -A 2 2 +l)/4 ^ _ ^ _ 2r) + ^ _ ^ + ^ ^ 

Because r is positive and A X A 2 > A only the second delta function contributes, giving the 
relation A = AiA 2 — It. As we are again considering the case when r < 1 our domain of 
integration for the other two variables is given by 1 < Ai < 1 + 2r and 1 < A 2 < ±±^i. When 
we perform the integral over A we are left with 

4r 2 (l - Ai 2 A 2 2 + 4rAiA 2 - 4r 2 ) 



2 A 2 2 - A x 2 - A 2 2 - 4r 2 ) 2 

x e -7r 2 x 2 (l+A 1 2 A 2 2 -A 1 2 -A 2 2 +4rA 1 A 2 -4r 2 )/4 _ ^ 

In order to evaluate this integral as a series in r it is useful to remove the r dependence from the 
limits. This is done by changing the integration variables using Ai = 1+ryi and AiA 2 = 1 + Ty 2 . 
Then the expansion of the parametric form factor is obtained by expanding the integrand for 
small values of r. 



K g ™(t,x) = I d Vl I dy 2 { - 2 V2 t + 



K 2(l- yi y 2 + yi 2 ) 2 
(2-y 2 )(2 + 2n 2 x 2 -6 yi -y 2 ) 
4(1 - 2/12/2 + J/i 2 ) 2 



y 1 (y 2 -2)(4-y 1 y 2 + 2y 1 2 -y 2 2 ) 
2(1 -2/12/2 + 2/1 2 ) 3 

t 2 + ...\. (16) 



Using Maple we performed this expansion up to seventh order and evaluated the integrals with 
the following result 

K goe (t, x) = 2t- (2ttV + 2)r 2 + (vrV + 2)r 3 - - ttV + ^ r 4 

/ttV 2ir 6 x e 2ttV \ 5 

+ 1 h 4 r 5 

12 3 3 J 

n w x w n 8 x 8 , , 7r 4 x 4 32 A , 

h ir 6 x 6 H 1 r 6 

60 4 3 5 J 

fn 12 x 12 n w x w 7n 8 x s 2n 6 x 6 ttV 32\ 7 

+ (^6F-^5- + ^2--^5- + ^ + yJ r ■ (17) 

For comparison with the semiclassical result it is convenient to extract an exponential factor 
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from this expansion. 



K GOE { 



T,X) 



+ 



-Bt 



5B 2 



2r - 2r l - (2B - 2)t a + [2B--)t 



3 

/4LB 3 



8B , - 

h 4 r 5 

3 1 



5B 2 



11B 2 



V 45 5 
where B has been defined after equation (11). 



+ 



32B 32 



AB + ?f) r 6 



r' + ... 



(18) 



4 Semiclassical approximation 

In this section we derive a semiclassical expression for the parametric spectral form factor. 
It is closely related to semiclassical approximations for the parametric two-point correlation 
function of the density of states [17, 18]. We start by expressing the density of states with the 
Gutzwiller trace formula [23] 

d osc (E,X) w -^-ReV A^exp ( ^S^] , where A, = ; 5 =e -^ /2 . (19) 

t» Y ^ / i? 7V /|det(M 7 -l)| 

The sum runs over all periodic orbits of the system with period T 7 , repetition number i? 7 , 
stability matrix M 7 and Maslov index /i 7 . 

The action is expanded in first order in the energy difference and the parameter difference 

( E ± 5 ± To x ± h) *> S ^ X) ± T ^ E - X) h ± Q -' {E ' x) To' (20 > 



where 



Q-, - ^ 



OX 



s = p M + ai (21 > 



is the parametric velocity. 

After inserting (19) and (20) into (6) and (8) and evaluating the integral to leading semi- 
classical order, we arrive at 

K(r, x) = ^r{J2 A,A;,e^^e^^5 [t - h±± J J , (22) 

where T# = 2nhd(E) is the Heisenberg time, and r = T/Th- Terms which have a sum of the 
actions in the exponent have been neglected, because they average away. The semiclassical 
expression (22) is the quantity that we will evaluate in the following. 

The diagonal approximation involves pairs of orbits that are either identical or related by 
time reversal. It has the form 

K(r,x) = ^(^\A^e^5(T-T^ , (23) 
where k is 2 if the system has time reversal symmetry and 1 if it does not. 
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One main ingredient in the following semiclassical calculation is the distribution of the 
parametric velocities Q 7 in the limit of very long periodic orbits. It has been shown that the 
Q 1 have a mean value of zero and a variance proportional to their period [16] 

(Q 7 )=0, (Q^-aT, T^oo, (24) 

where the averages are performed over trajectories with period around T. It has been motivated 
that the Q 1 have a Gaussian distribution for long periodic orbits (see [18]) and this is the main 
assumption that we will use in the following. The proportionality factor a in (24) is related to 
the variance of the level velocities. [24, 21] 

° 2 ~ a 4 ■ < 25 > 

Using the Gaussian assumption for the distribution of the parametric velocities and (25) we 
perform the average over the Q 7 , assuming that it can be done independently, and obtain 

( e — ) = e"^^ = e ~ BT/TH , (26) 

where B = 2tt 2 x 2 /k, is the same quantity as defined after equation (11). The remaining sum 
over periodic orbits can be evaluated with the Hannay-Ozorio de Almeida sum rule [5] 

|A 7 | 2 <5(T — T 7 ) ps T . (27) 

7 

We find that the diagonal approximation is given by 

K(t) = nre- BT (28) 
in agreement with the first term in the expansion of the random matrix results, (12) and (18). 



5 Off-diagonal contributions 

The off-diagonal terms of the parametric form factor are due to pairs of trajectories which are 
correlated [7, 10]. In the following we briefly review the main steps in the derivation of the 
semiclassical expansion of the spectral form factor (in our notation K(t,x = 0)) according to 
[10, 11]. The correlations that are important for the expansion of the form factor for small r 
are due to close self-encounters of a periodic orbit in which two or more stretches of an orbit 
are almost identical, possibly up to time reversal. In general, a long periodic orbit has many 
of these encounter regions, and they are connected by long parts of the orbit, the so-called 
"loops". The correlated pairs of orbits are almost identical along the loops, but they differ 
in the way in which the loops are connected in the encounter regions. Correlated orbit pairs 
have certain "structures" that are characterized by the number of encounter regions V in which 
the loops are connected in a different way, the number of involved orbit stretches l a in each 
encounter region a, and the way in which the loops are connected by these stretches. A more 
accurate definition of structures can be given by putting them in a one-to-one relation with 
permutation matrices that describe the reconnections of the loops. One defines further a vector 
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v whose l-th component, vi, specifies the number of encounter regions with / stretches, and the 
total number of orbit stretches is denoted by L. Hence 

V = 5>,, L = ^Z a = X>. (29) 

l>2 a l>2 

The semiclassical contribution to the form factor is evaluated in two steps. First the summation 
over orbit pairs with the same structure is evaluated by using that long periodic orbits are 
uniformly distributed over the surface of constant energy in phase space. Then the summation 
over the different structures is performed which is a combinatorial problem. 

In the following we present some details of this calculation. In each encounter region a with 
l a orbit stretches one chooses a perpendicular Poincare surface that is centered on one of the 
stretches. The relative positions of the piercings of the other stretches through the Poincare 
surface are described by coordinates along the stable and unstable manifolds. The partner 
periodic orbit connects the loops that start and end at the encounter region in a different way 
and the resulting contribution to the action difference is given in the linearized approximation 
by 

i a -i 

(AS) a = ^ S *jUaj , (30) 
3=1 

where s a j, u a j, j — 1, . . . , l a — 1 are appropriate differences of the coordinates along the stable 
and unstable manifolds. For ease of notation we discuss here the two-dimensional case in which 
the coordinates s a j and u a j are scalars. If the Poincare surface is moved along the stretches 
in the encounter region these coordinates decrease or increase, exponentially, however their 
product remains constant. The length of the encounter region is determined by requiring that 
all coordinates remain smaller than an arbitrary small constant c whose exact size is not relevant 
for the following calculations. 

The uniform distribution of the long periodic orbits on the energy shell is then invoked to 
sum over all orbit pairs with the same structure. It is convenient to also sum over all structures 
with the same vector v and express it in the form 

fixed v r / \ 

K v (t ) = ±- Y \A,\ 2 e iAs ^ h 5(T - T 7 ) = N{v)kt / d L ~ v s d L ~ v u Wt V> u > J™/* (31) 
Th . , . J L 

(7,7 ) 

so that K(t) = kt + K v (t). N(v) is the number of structures with the same v, and s and 
u are vectors whose components are the s a j and u a j for all a and j. In (31) the amplitudes 
and periods of the two correlated orbits are set equal. wt(s, u) is the density of the self- 
encounters for a given structure and separation coordinates s a j and u a j. For long orbits it is 
given asymptotically by 

w T (s,u) _T(T-E a U° ' 1 



one / 



L Ll^-VUA " (32) 

The factor 1/L in (31) takes care of an overcounting related to the choice of an initial point of 
the trajectory [10]. 

The integral in (31) is evaluated by using 



J n dsa J 




d Uaj « { ^ - -i<**>i_ (33) 
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Hence after expanding the numerator of (32) the only term that survives is the one that contains 
a product of all encounter times t" nc which is cancelled by the denominator. Therefore one can 
replace the density wt in the integral by 

L L\VL L ~ V (L-V-l)\ ' 1 ' 

where we used that Yla^a = Yli^ 1 - Inserting (34) into (31) and evaluating the integral with 
the use of (33) yields the contribution from orbit pairs with fixed v, 

KM = KT .-v +lN{v) (-WW^ . (35) 

As can be seen, vectors v with the same value of L — V + 1 contribute to the same power of 
r. Hence the power series expansion of the form factor can be obtained by summing over all v 
with the same L — V + 1. 

oo 

K(r) =KT + J2 K nT n , (36) 



n=2 



where 



L~V+l=n 

K n = 7 -^— N(v) = (N(v)(-l) v Hn/L. (37) 

V I 



71 



By deriving a recurrence relation for the numbers N(v) it is shown that the semiclassical 
expansion in (36) agrees with the small r expansion of the form factor of RMT in cases with 
or without time reversal symmetry [9, 10, 11]. 

Let us now come back to the parametric form factor K(t, x). We have to evaluate 

K(r..v) = — (> \ \ ^' ) c^S(:r--lv\ . (38) 




It now contains an additional term involving the parametric velocities Q 7 - As for the diagonal 
approximation we assume that the average over the Q 1 can be performed independently from 
the actions and amplitudes of the orbits. However, we have to look at the average (26) more 
carefully. It is valid as long as there are no systematic correlations between different parts of 
the same periodic orbit. Long periodic orbits have many encounter regions in which different 
stretches of an orbit are almost identical and cannot be considered as independent or uncor- 
rected. When we calculate the contribution to the form factor of orbit pairs with a certain 
structure, then we average over periodic orbits that all have the same number and types of 
encounter regions. The corresponding correlations between different parts of an orbit have to 
be taken into account. We should consider the average over the parametric velocities for the 
loops and encounter regions separately. The contribution from a loop is 



ixQ 



Loop 



<e-i*-> = e ~ BT ^ /TH , (39) 
while the contribution from the / orbit stretches in an /-encounter region is 



e ^- } = e -BiH mc /T H (4Q) 
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This means that the average over the parametric velocities is now given by 

( e ^) = e -B(T-Z a l a t^ c )/T He -BZ a l 2 a t^ c /T H _ ^ 

The contribution of all orbit pairs with encounter regions described by v then have the form 



fixed v 



Kv ^ x) = J_ V \A,\ 2 e iAs ^ h e^5(T - T 7 ) = iV( V ) K r / d £ - v u £li£l^ e is ^ , 

-iff ^— ' y l 

(7:7 ) 

(42) 

where 

zt(s,u) w t (s,u) ^ e- OT ^T(r - E Q UL)^ 1 EL e-*^- 1 ^* 

Here a labels the V different encounters, each being a /^-encounter. Again because of (33) the 
only terms that contribute in the semiclassical limit are those where the encounter times in the 
numerator and denominator cancel exactly. As a first step we can expand the exponentials as 
a power series up to first order 

z T (s, u) _ e~ BT / T »T(T - £ a U^" 1 EL(1 - I«MBM _ {u) 



un L - v nt 

1 let 



a 

enc 



To obtain a product of the V different encounter times in the numerator we can take r of them 
from the product over a and V — r of them from the bracket with the exponent L — 1. The 
corresponding coefficient is obtained by combinatorial considerations. Then we sum over all 
values of r from to V, and the result is 

Z T (S,U) e -BT/T HT V T L-V + r-l B r {L _ m _ 1)V ft g. 



(45) 



L L\n L ~ v ^ 7^(L-y-l + r)!H 

r=0 " v y oi,...,a r 

distinct 

We insert this into (42), evaluate the integral with formula (33), and obtain 

V ~r /• • oi,...,o r 

distinct 

The contribution of orbits for the different types of encounters can be calculated with this 
formula, and they are shown for orbit pairs with L — V < 4 in table 1. The vectors v are 
represented in the form (2) V2 (3) V3 . . . and the horizontal lines separate vectors v with different 
value of L — V. The numbers N(v) can be calculated by combinatorial methods [11]. 

To find the total contribution to the form factor we now multiply the middle column that 
contains K v /(kN(v)) by k and N(v), add the diagonal approximation and sum over different 
v. If we do that for all orbits pairs with L — V < 8 for the case without time reversal symmetry 
(k — 1), we obtain the expansion for the form factor in r up to 9th order 



B 2 r 5 /iV 
K(T)=e-» T + — + — + ... 



(47) 
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V 


L 


V 


K v (t,x)/(kN(v)) 


7V(v), no TRS 


N(v), TRS 


(2) 1 


2 


1 


-e~ Br (r 2 + .Br 3 ) 




1 




A 


o 
_ 


p -Bt ( 3 1 d_4 , S 2 r s \ 


1 
1 





(3) 1 


3 


1 


_ e -Br ( T 3 + 5r 4) 


1 


4 


(2)3 


6 


3 


Bt / 2r 4 | 2Br 5 . B 2 r 6 , B :i r 7 \ 
C ^ 3 ~ l ~ 3 1 6 1 90 y 


- 


41 


(2) 1 (3) 1 


5 


2 


n -Br ( 3r 4 , 3Bt 5 , B 2 r 6 ^\ 

^ 5 5 1 10 y 


_ 


60 


(4) 1 


4 


1 


-e- Br (^ + ¥) 
v z z / 


_ 


20 


(2) 4 


8 


4 


Bt ( r 5 , Br 6 , B 2 r 7 , B 3 r 8 , B 4 r 9 \ 
I 3 3 10 1 90 1 2520 ) 


21 


509 


(2) 2 (3) 1 


7 


3 


Bt ( 2r 5 , 2Br 6 , B 2 t 7 , B 3 r 8 \ 
e 17^7 1 14 1 210 / 


49 


1092 


(2)W 


6 


2 


„-Bt / 2r 5 | 2Br 6 , B 2 t 7 \ 
C ^9^9 1 30 J 


24 


504 


(3) 2 


6 


2 


p -Bt (t 5 , Bt« B 2 T 7\ 
e ^ 4 T 4 T 20 j 


12 


228 


(5) 1 


5 


1 


-<*-* a + ¥) 


8 


148 



Table 1: Contribution of different types of orbit pairs to the parametric form factor 



This agrees with the first three terms of the expansion (12) in the section on RMT. It is 
noticable that when summing over terms with the same value of L — V, that all terms cancel 
apart from the highest order term from orbit pairs with only 2-encounters. In fact we will show 
this using a recurrence relation in appendix A. This allows us to calculate the expansion of the 
form factor to all orders in r which is done in the next section. 

For systems with time reversal symmetry (k = 2) we sum over all contributions with L — V < 
6, and obtain the expansion of the parametric form factor in r up to 7th order 



K(t) = e 



-Bt 



2t 



2r 2 - {2B - 2)r 3 + ( 2B - ) r 4 



+ 



5£? 2 



3 

415 3 



T 
11s 2 



- !|r 5 - 
325 



+ 



5B 2 n 32 
4B + — 

3 5 



32 



45 5 5 

This agrees with the expansion of the RMT result (18). 



T lr' + ... 



(48) 



6 Systems without time reversal symmetry 

In this section we derive the full expansion of the parametric form factor for small r for the 
case of systems without time reversal symmetry. For this purpose we rewrite the expansion 
K(t,x) = re~ Br + Yl°^=2 Kv{t, x) with K v (t,x) given in equation (46) in the following form 

00 n— 1 

K(r, x) = re~ BT + e~ Br £ £ K n ,,r n+r B r . (49) 

n=2 r=0 
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The coefficients K nr have the form 

Kn,r = J^y S n[fr(v)\ , (50) 

where 

L-V+l=n 

Sn[fr(v)]= fr(v)N(v), (51) 

v 

and the functions f r (v) are given by 

frW= a iV-V+rVH ^ (i ai -l)x...x (/«,-!). (52) 

V 'J cei,...,a r 

distinct 

The first few functions are 
fo(v) = 1 , 

f M _ E«(*« - 1) _ EkMk-i) _ 

Jl{ ] L-V L-V ' 

f / \ _ E ai , a2 (^-i)(^-i)-EJ^-i) 2 (L-\/) 2 -E fc ^(fe-i) 2 , , 
/2W 2(L-y)(L-y + i) 2(L-y)(L-y + i) ' 1 j 

We need to evaluate the quantities S n [f r (v)] for r < n. In appendix A it is shown that 
Sn[fr(v)] — for r < n — 1. Hence the only non- vanishing terms in the expansion (49) are 
those with r — n — 1. 



oo 



X(r, x) = re- ST + ^ T^r 2 "- 1 ^ 1 . (54) 

n=2 

Since r satisfies r < V we have V > n — 1. Together with the condition L - V = n - 1 we find 
that 2V > L. This is only satisfied for orbit pairs with V 2-encounters for which v = (2) y and 
L = 2V. The contribution of these orbit pairs to the form factors can be calculated explicitly. 
We obtain from equations (50), (51) and (52) with r — V — n — 1, L — 2V, and l a — 2 for all 

a, 

N(v) (L-V-l)\ (-l)"-^"- 1 
= (n-2)!(L-V-l + r)! = (2n - 2)! ^ " (55) 

The number JV(u) can be obtained from an explicit formula that has been derived for systems 
without time reversal symmetry in [25]. In our notation it has the form 



The notation i/ < v means that the sum runs over all integer vectors v' whose components 
satisfy < v ' k < v & for all k. Furthermore, V = L(v') and V = V(v'). In the case of vectors 
v of the form (2) v the only non- vanishing component of v is v 2 = V and the sum runs over all 
vectors with component v' 2 = m where m — 0, . . . , V. The result is 

1 ^ (2m)\ (2V - 2m)\ (2V)\ (-1)^ + 1 

{) 2V + 1^ 2 v m\(V-m)\ 2 V (V + 1)\ 2 ' 1 ' 

m=0 v ' v ' 
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where the last equality can be found in [26] . The expression vanishes if V is odd or equivalently 
if n is even. With V — n — 1, n odd we obtain K n n _\ = 1/n! and complete expansion of the 
form factor is (n — 2k + 1) 

in agreement with the RMT result (11). 



7 Conclusions 

This work is a continuation of recent developments in semiclassical periodic orbit expansions. 
These methods have been applied to several spectral and transport quantities in order to 
demonstrate the universality of quantum fluctuation statistics of chaotic systems. We extended 
these ideas to obtain a semiclassical expansion of the parametric spectral form factor K(r, x) 
for small r in agreement with RMT. For the GUE case we showed agreement for all orders, 
while for the GOE case we showed agreement for the first seven terms in the expansion. These 
terms can actually be obtained more quickly by the semiclassical method than by an expansion 
of the double integral for the GOE result. 

The main input that is needed in addition to the semiclassical calculation of the spectral 
form factor K(r) is the distribution of the parametric velocities of long orbits. This is the 
distribution of the derivatives of the actions of periodic orbits with respect to the parameter, 
which is commonly assumed to be Gaussian. In addition, when averaging over orbit pairs whose 
correlations are described by a particular structure, one has to take into account correlations 
between different parts of the same trajectory. These correlations are due to almost identical 
orbit stretches within the encounter regions. 

The limitations of the semiclassical calculation are similar to that for the spectral form 
factor. One main open point concerns the region r > 1. In this regime the random matrix 
expressions for K(t, x) have a different functional form. So far it is not known how to extend 
the semiclassical approach to this regime except for a small region near r — 1 [17, 18]. 



A Recurrence relations 

In this appendix we show that the quantities S n [f r (v)], defined in (51) and (52), vanish for 
r < n — 1. The function f r (v) are defined in terms of a restricted sum in which all summation 
indices are distinct. As a first step this sum is expressed by unrestricted sums. How to do this 
by a combinatorial sieving is discussed, for example, in section 4 of reference [27]. 

We first introduce some notation. A set partition F of the set of integers {1,2,... r} is a 
decomposition of this set into disjoint subsets [F 1: . . . , F v \. Then \F-i\ + . . . + \F V \ = r where 
\Fi\ is the number of elements in the set Fj. Let us define a generalization of the Kronecker 
delta-function 

F J 1 if (Xi = aj for all % and j such that i,j G Ft for some k, 

d «w* = i otherwise. ( } 
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Then 

E [■■■] = E^) E CJ-i- (so) 



ai,...,a r jr ai,...,a r 

distinct 



where the first sum of the right-hand side runs over all set partitions of the set of r integers, 
and the corresponding Mobius function is given by 

M^) = f[(-l) |Flhl (l^|-l)!- (61) 
i=i 

If we apply this to the functions f r {v) we obtain 

MV) = (L -V-l + rjl r! £ "^'M ' (62) 

where 

9F(v)^(^v k (k-ir^x...x^v k (k-ir^ . (63) 

The first few functions f r (v) were given in equation (53). 

The expansion of the form factor K(t) was evaluated in [9, 10] by using recurrence relations 
for the number of structures M{v) corresponding to a vector v. These recurrence relations were 
obtained by relating orbits with L loops to orbits with L — 1 loops by considering all possible 
ways of removing a loop (i.e. letting its size shrink to zero). 

For systems without time reversal symmetry the relevant recurrence relation is 

v 2 N(v) + J2 v k k ^ k+1] kN(v^ k+ V) = . (64) 

k>2 

Here the vector v \ k ^ k + l \ [ s obtained from the vector v by decreasing the components v k and v 2 
by one and increasing the component v k +i by one. Hence L(v^- k ' 2 ~* k+1 ^) = L(v) — k — 2 + (k+l) = 
L(v) - 1 and v(v^ 2 ~* k+ ^) = V(v) - 1. 

In order to obtain the coefficient of the form factor expansion one has to sum over the 
numbers Af(v) for all vectors for which L(v) — V(v) + 1 = n. The recurrence relation may be 
used for this purpose, because one can show that for each k 

L-V+l=n L'-V'+l=n 

v{%- k+1] h(v [k ^ k+1] ) = E < + Mv'), (65) 

v v' 

where h(v) is some function of v. One condition is that v\ = v ^ k ' 2 ^ k+1 ^ = o, because the vectors 
describe encounter regions which contain at least two orbit stretches. Summing the recurrence 
relation (64) over v yields 

= S n [v 2 + v k+ik] = S n [L -V] = (n- l)S n [l] . (66) 

k>2 

This shows, for example, that all off-diagonal terms of the form factor K{r, 0) vanish [9, 10]. 
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We want to show in the following that S n [gF(v)] — if r < n — 1. We consider first the case 
when the partition consists of only one subset F 1 with \Fi\ = r. Then gp(v) = g r (v) where 

g r (v) = J2Mk-l) r ■ (67) 

k 

We show that S n [g r (v)\ = 0ifr<n — lby induction. The statement is true for r = 0, because 
S n [l\ = by equation (66). Now we fix a value of r < n — 1 and assume that the statement is 
true for all smaller values of r. From the definition (67) follows that 

g r (v [k ^ k+1] ) = g r (v) -h r (k) , h r (k) = (k-l) r -k r + 1 . (68) 

Points that will be important in the following are that h r (l) = and that h r (k) is given by a 
finite power series in k whose highest order term is —rk r ~ l . 

Multiplying equation (64) by g r (v) and using relation (68) we obtain 

= v 2 g r (v) N(v) + J2< + ik9r(v')N(v') +Y,v' k+1 kh r (k)N(v') , (69) 

k>2 k>2 

where v' = v i k ' 2 ^ k + 1 }, i n the last sum we can start the sum at k — 1, because h r (l) = 0, and 
then change the summation index k — > k — 1. After summing over all vectors v we obtain 

= S n [v 2 g r (v) + ^v k+1 kg r (v) + ^v k {k - l)h r (k - 1) 

k>2 k>2 

= S n [(L - V)g r (v) - J2^r(k - If + . . .] . (70) 

k>2 

In the second line we used that v 2 + X]fe>2 v k+\k = X^>2 v i(h ~ = L — V, and we wrote only 
the highest order term of h r (k— 1). The lower order terms, denoted by the dots, involve powers 
(k — l) m with m < r and can be neglected due to our induction assumption. Hence we find 
that 

(n-r-\)S n \g r (v)\=Q, (71) 

so that indeed S n [g r (v)] — if r < n — 1. The proof for general gp(v) is very similar. We 
consider the general form 

V 

9F(v) = l[g\ Ft \(v) , (72) 

i=l 

and we use again induction to prove that S n [gF] — if r < n — 1. The statement is true for 
r = 0, and we fix a value of r and assume that it is true for all smaller values of r. In order to 
use the recurrence relation (64) we note that 

gAv)^f[(9\F 4 (v^ k ^) + h lFil (k)). (73) 
i=\ 

We multiply equation (64) by gp{v) and use relation (73) to obtain 
=v 2 g F (v) N(v) + v' k+1 kg F (v')N(v') 

k>2 

v v 

+ + V#)) - l[9w(v')}N(v') , (74) 

k>2 i=l i=l 
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where we added an additional term and subtracted it again. As before v' = ^[ fe > 2 - >fe + 1 ]. In the 
second sum we can start the sum ctt h — 1, because hi(l) = for all %, and then change the 
summation index k — > k — 1. After summing over all vectors v we obtain 

V V 

= S n [v 2 g F (v) + ^2v k+1 kg F (v) + ^2v k (k - l)[Y[(g\ Fi \(v) + h m {k - 1)) - 

k>2 k>2 i=l i=l 

= S n [(L - V)g F (v) + J>(* - 1) Y,{~\F 3 \{k - l)^'" 1 ) \{g m {v) + ...]. (75) 

k>2 j=l i^j 

In the step from the first to the second line we expanded the first product, inserted the power 
series for the functions h\ Fi \{k — 1) and wrote only those terms that do not vanish due to the 
induction assumption. We obtain further 

V 

= S n [(L - V)g F {v) - \ F i\ m\(v) l[g\ Fi \(v) + ...] 

= (n-l-r)S n [g F (v)}, (76) 
which concludes the proof that S n [g F (v)] = for r < n — 1. 
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